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Abstract 



The dependence of fusion cross section on the isotopic composition of col- 
hding nuclei is analysed within the dinuclear system concept for compound 
nucleus formation. Probabilities of fusion and surviving probabilities, in- 
gredients of the evaporation residue cross sections, depend decisively on the 
neutron numbers of the dinuclear system. Evaporation residue cross sections 
for the production of actinides and superheavy nuclei, listed in table form, 
are discussed and compared with existing experimental data. Neutron-rich 
radioactive projectiles are shown to lead to similar fusion cross sections as 
stable projectiles. 
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Typeset using REVTgX 



The synthesis of superheavy elements (^=106-112) was reached by cold fusion of heavy 
ions with lead and bismuth targets Hot fusion reactions using ^^^Th, ^ss^j g^j^^ 242, 244 
targets were also applied to synthesize the elements with Z=110, 112 and 114 0. A pos- 
sible next step is to explore the synthesis of heaviest nuclei with radioactive beams 
Microscopical investigations of such planned experiments are a challenge for theory. Usu- 
ally the surviving probability Wgur of the formed compound nucleus against fission in the 
de-excitation process is considered as the crucial factor which is mainly responsible for pro- 
ducing heavy and superheavy elements. With neutron-reach projectiles one can obtain a 
large stability (large Wgur) of the compound nucleus. However, the probability of complete 
fusion PcN, dependent on nuclear structure effects and on the neutron excess above the 
nearest closed shells in the colliding nuclei, is also very important for the correct calculation 
of the evaporation residue cross section aER- For example, experimentally extracted proba- 
bilities PcN are strongly decreased when the neutron numbers of the projectile or target 
deviate from magic numbers. 

The existing fusion models can be distinguished by their choice of the relevant collective 
degree of freedom responsible for complete fusion. Many models assume an adiabatic melting 
of the nuclei along the relative distance R of nuclear centers (or the elongation of the 
system) |[7|-p!^. However, it was demonstrated that the adiabatic scenario of fusion along 



the relative distance leads to a large overestimation and an incorrect isotopic trend of the 



fusion probability [T^. The dinuclear system (DNS) concept |Tj-|T9[ assumes that the united 
system is reached by a series of transfers of nucleons or small clusters from the light nucleus 
to the heavier one in a touching configuration. So, the dynamics of fusion is considered as a 
diffusion of the DNS in the mass asymmetry, defined by 77 = {Ai — A2)/{Ai + A2) {Ai and 
A2 are the mass numbers the DNS nuclei), where the potential barrier B^^^ in r] supplies a 
hindrance for fusion. 

The assumption of a touching configuration of the two reacting nuclei in the DNS model 
is supported by the structural forbiddenness of fusion p(]| , pT| which hinders the nuclei to 
melt together along the relative distance. This aspect is phenomenologically described with 
a double folding potential in frozen density approximation which shows a minimum near 



the touching distance of the nuclei There are also experimental evidences p3| , p^ that 



the mass asymmetry degree of freedom equilibrates more rapidly than the elongation of the 
system. During the characteristic time of fusion a statistical approach is applicable to treat 
the evolution of the DNS which also includes a diffusion to larger relative distances between 
the centers of the nuclei describing the quasi-fission process (decay of the DNS) competing 
with the complete fusion. In reactions with heavy nuclei the quasi-fission channel dominates 
and leads to a strong reduction of the fusion []T3|-[T^. 

In accordance with the DNS concept the evaporation residue cross section is factorized 
as follows ||17| 



(TER{Ecm) = (rciEcm)PcNiEcm, J = 0)Wsur{Ecm, J = 0). (l) 

The calculations of the evaporation residue cross sections demands an analysis of all three 
factors in ([^). The value of cTc is the effective capture cross section for the transition of the 
colliding nuclei over the entrance (Coulomb) barrier with the probability T: 

a,{E,^) = TlX^Jma. + l)'T(Ee^, J = 0). 
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The contributing angular momenta in the evaporation residue cross section are limited by the 
surviving probability Wsur{Ecm, J) with Jmax ~ 10 — 20 when highly fissile superheavy nuclei 
are produced This corresponds to almost central collisions with impact parameters 

smaller than 1 fm. The value of Jmax is smaller than the critical angular momentum Jcrit 
which restricts the capture cross section. For reactions leading to superheavy nuclei, values 
of Jmax = 10 and T{Ecm, J = 0) = 0.5 are chosen for energies Ecm near the Coulomb barrier. 
The capture cross sections obtained with these parameters are in agreement with the ones 
calculated within a microscopical model [p!^] . 

The probability of complete fusion Pcn in (0) depends on the competition between 
complete fusion and quasi-fission after the capture stage. It can be expressed by rates in 
the quasi- stationary follows 
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As in Ref. [0 we use a two-dimensional Kramers-type expression (quasi-stationary solution 
of the Fokker-Planck equation) with the quasi-stationary rates of fusion X^^ and quasi-fission 
through the fusion barrier (-Bfug) in rj and quasi-fission barrier (-Bqf) in R, respectively. 
The second term in (0) is related to the transient times r^j and to reach the quasi- 
stationary rates along the R and 7] coordinates (r^^, > A;^'', A^'') |16| . 

In the case that the fusion barrier is much higher than the quasi-fission barrier, B^^^ ^ 
-Bqf, i.e. if the transient time in t] is larger (or equal) than the lifetime to of the initial 
DNS, we obtain IITBI 



CN 



\Kr 
1.72 



r^(exp[to/r^ 



1) -tn 



(3) 



Since the pocket in the nucleus-nucleus potential becomes very shallow (-Bqf ~ 0) in reac- 
tions with large Zi x Z2, the lifetime to of the DNS is strongly depressed with increasing 
bombarding energy Ecm. above the Coulomb barrier. Due to this, the value of Pcn in Eq-(3) 
is smaller than the one in the quasi-stationary regime, given by Eq. (2), which can not be 
reached in this case |]16| . 

The surviving probability under the evaporation of x neutrons is considered according 
to [0,0 as 



WsuriE*cN,J)-PUE*CN,J)Il 



^niE^^ , Ji^ 



1 ^n{EcN,-, Ji) + ^f{EcN,^ Ji, 



(4) 



Here, P^n is the probability of realization of the xn channel at the excitation energy E'^^ 
of the compound nucleus, i the index of evaporation step, r„ and V f are the partial widths 
of neutron emission and fission. -E^TVi ^"^^ Ji mean values of excitation energy 

and angular momentum quantum number of the compound nucleus, respectively, at the 
beginning of step i with Eq^^ = -E^Af ^.nd Ji = J. In the calculation of Wsur we used 



the microscopic corrections of MoUer and Nix |2^ as fission barriers. The neutron binding 



energies are also taken from [27 



The barriers -B^^g and -Bqf are given by the potential energy of the DNS which is calculated 
as the sum of binding energies -Bj of the nuclei (i=l,2) and of the nucleus-nucleus potential 
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U{R,v, J) = Br + B2 + V{R,v, J) - [5i2 + vUJ)]- (5) 

The shell effects are included in the binding energies. The isotopic composition of the 
nuclei forming the DNS is obtained with the condition of a A^/Z-equilibrium in the system. 
The value of U{R, i], J) is normalized to the energy of the rotating compound nucleus by 
Bi2 + Vrof Deformation effects are taken into account in the calculation of the potential 



energy surface |17]. The heavy nuclei in the DNS, which are deformed in the ground state, 



are treated with the parameters of deformation taken from Ref. The light nuclei of the 
DNS are assumed to be deformed only if the energy of their 2+ state is smaller than 1.5 
MeV. As known from experiments on sub-barrier fusion of lighter nuclei, these states are 
easily populated. For the collision energies considered here, the relative orientation of the 
nuclei in the DNS follows the minimum of the potential energy during the evolution in 77. 

The experimentally observed hindrance of the fusion roughly increases with growing 
Coulomb repulsion between the colliding nuclei, but also their shell structure and isotopic 
composition play a major role P,|29|-pT|. In Table 1 we present calculated excess energies 
above the entrance Coulomb barrier in the DNS model for various reactions and compare 
them with the surplus of energy extracted from experimental data above the corresponding 
Bass barriers [Q. In these calculations we did not average the inner fusion barrier B^^^ 



over all possible orientations of colliding nuclei as we usually do in the calculations of Pcn 
and (Ter, taking approximately the half of the deformation parameters of the nuclei of the 
entrance channel. The obtained energy thresholds in t] are maximal ones. They are not 
always in good agreement with the data extracted from experiment because these data are 
not directly measurable but are obtained with model assumptions about Pcn and Wgur- 

As shown in Table 1 the isotopic trends of the DNS model agree with the experimental 
ones. The energy thresholds for fusion increase and, correspondingly, the fusion probabilities 
decrease [ pU| , pn| when the neutron number of projectile or target deviates more from a magic 



number in the reactions 902r+90,92,962r, 90,962r+i00Mo 86Kr+99.i02,i04j^^^ 90,92,94,962^+1243^ 
and ^^Kr+^^'^'^^^Xe. This effect is simply explained by the deformation of the nuclei in the 
initial DNS and DNS at the top of the barrier in 77 and by the shell effects in dependence of 
the DNS potential energy on 77. For example, the value of the energy threshold for fusion, 
which determines the fusion probability, is larger in the ^^Kr+^^°Xe reaction than in the 
^^Kr+^'^^Xe reaction Since in addition the surviving probability Wgur is larger in the 



reaction with ^^^Xe than in the reaction with ^^°Xe, there results an experimental difference 



of about 3 orders of magnitude in aER in these reactions [^. For most reactions, for 
example, ^°Zr+^°Zr, ^°''Mo+^°''Mo and ^^'^Pd+^^^Pd, we obtained evaporation residue cross 
sections with the values of Pcn of the DNS model which are in good agreement with the 



experimental cross sections |T5|J16f| . In contrast, models which treat fusion as a motion in 
R, give an incorrect isotopic trend of Pcn- In these models Pcn always increases with the 
neutron number above the nearest closed shell |]6|,|32[| because an increasing deformation of 



the nuclei in the entrance channel effectively lowers the barrier. 

Since the evaporation residue cross section increases with the number of neutrons in all 
reactions listed in Table 1, the value of Wsur has to grow faster than Pcn decreases. In 
fusion reactions leading to actinides, for example in the ®^'™Zn+^^^Yb reaction, the increase 
of Wsur with the neutron number of the system is stronger than the decrease of Pcn This 
effect, shown in Fig. 1 for reactions ^Zn+^^^Yb, gives a certain preference for neutron-rich 
projectiles in producing actinides. Note that the numbers of neutrons in the nuclei ^^Zn and 
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^^Zn are close to different magic numbers. 

Fusion probabilities in symmetric and almost symmetric reactions with heavy nuclei like 
i24,i32gj^ and ^^^'^^^Xe strongly depend on the model of fusion. For example, in an adiabatic 
treatment, where the fusion is mainly described by the dynamics in the relative distance 
coordinate with an increasing neck, we found Pcn ~ 10~^ and aER ~30pb for the reaction 
^^^Sn + ^^^Sn— *-^^^Fm+3n. In the DNS model the values of Pcn and, correspondingly, aER 
are about three orders of magnitude smaller. According to the DNS concept, cross sections 
for the synthesis of the heaviest elements in nearly symmetric reactions are very small due 
to small fusion probabilities, for example, in the reactions ^^^Xe + ^^^Xe— >-^^^Hs, ^^^Xe + 
150,154^^^292,296^14^ i32gj^ ^ 160(^^^292^14 137^6 + ^^^^ui-^'^'^^'lU. Experimental data 

on symmetric reactions with stable and radioactive beams could help to prove the DNS 
model for the fusion process and would give information about the time for the transition 
from the diabatic to adiabatic regime (the time of suppression of the structure forbiddenness 



for melting of nuclei |20]). 

In contrast to other models, the optimal excitation energy -E^Af compound nucleus 

and evaporation residue cross section aER in cold fusion reactions with stable projectiles 



are reproduced in the DNS concept |jT^. These results are listed in Table 2 for reactions 



leading to the Fm element and Pb- and Bi-based reactions. The evaporation residue cross 
sections are compared with the experimental data of Refs. [0,^]. All other cross sections 
are predictions of the present version of the DNS model. Figs. 2a) and 2b) show the fusion 
probabilities and the optimal excitation energies of the compound nuclei, respectively, for 
208pb, 209Bi(^X, In) reactions. 

The values of the optimal excitation energy -E^Af calculated by applying theoretical Q- 
values of Refs. |27]. They increase for Z > 112. Q- values of Ref. [|4] are slightly different for 



Z > 113. As in the case of reactions with heavy nuclei mentioned above, the calculated values 
of Pcn are maximal when the neutron number of the projectile is equal to a magic number, 
for example, in the reactions ^^Ge+^^^Pb, '^^Se+^^^Pb and ^^Kr+^^^Pb. The decrease of the 
cold fusion cross section by four orders of magnitude from Z =104 to 112 is mainly caused by 
a decrease of Pcn due to a strong competition between complete fusion and quasi-fission in 
the DNS (see Fig. 2a)). For the reaction ^°Zn+^°^Pb^^'^'^112+ln, aER ~lpb is practically on 
the level of the present experimental possibilities. In reactions '^^''^^Ge+'^^^Ph-^'^^^''^^^114:+ln 
we expect a value of aER which is smaller than 0.2 pb. The values of aER for the Z=11Q and 
118 elements formed in the ^^Se, ^^Kr+^^^Pb reactions are of the order of 0.01 pb (Table 2). 

The values of Wsur in Table 2 were calculated with the theoretical data of Ref. |]2^. One 
can see that characteristic values of Wsur are about 10^^ — 10"^ for nuclei with Z=104-113 
and about 10~^ for nuclei with Z=114, 116 and 118. The proton magic number 114 in the 



region of the stability island [ 27, 34, 35 1 leads to a larger increase of Wsur- The large Ws- 



of the nuclei ^^^114, ^^^116 and ^^^xxg arises due to the fact that the neutron number in 



these nuclei is equal to the theoretically predicted magic number A^=178 [p7| , p^ . When 
the number of neutrons deviates from this magic number, Wsur decreases. The surviving 
probabilities in the reactions ^''Zn+^'^^Pb and ^^Ge+^°^Pb were calculated with the data of 



Ref. 0] because Wsur becomes unrealistically small with the data of Ref. ||2^ (about two 



orders of magnitude in comparison to the neighbouring nuclei). Using the microscopical 
corrections and neutron binding energies from Ref. ||3^ instead of Ref. we obtained 
even smaller cross sections aER of 50, 5 and 0.8 fb for the reactions with the projectiles 



5 



76 



Ge, ^^Se and ^^Kr on ^°^Pb, respectively. With the fission characteristics of Ref. [p5 



we 

get evaporation residue cross sections aER which are again smaller than the latter ones. In 
conclusion the cross sections in Table 2 are optimistic estimates. 

Let us consider whether the expected values of evaporation residue cross sections are 
larger with radioactive projectiles. In the Pb-based reactions with neutron-rich nuclei 
7o,74,78]\j-j^ ®°Zn, '''^~®^Ge, ^^~^^Se and ss-s^j^^, ^j^g inner fusion barrier B*^^^ in mass asym- 
metry varies between 12 and 22 MeV. In order to overcome this barrier, the initial DNS 
must have excitation energies which lead to an excited compound nucleus with the possi- 
bility of an In or 2n emission. The calculated cross sections for some possible reactions 
are presented in Table 2. Bombarding energies near the Coulomb barrier lead to maximal 
evaporation residue cross sections. In these reactions the increase of Wsut is compensated 
by a decreasing fusion probability Pcn and the value of aER depends weakly on the isotopic 
composition of the colliding nuclei. The values of Pcn and Wsur for the reactions "^Ni+^^^Pb 
and "^Ge+^°^Pb are presented in Fig. 3 as functions of A. The calculations were performed 
with the same parameters as used for the stable projectiles and are in good agreement with 
experimental data (Table 2) |I|,p!7|. Due to deformation effects and binding energies of the 
nuclei in the DNS, the dependence of Pcn on ^ can have some minima and maxima. 

The yield of the element Z = 110 results larger in the ^^Ni+^°^Pb reaction than in 
the 62,64]\jj_|_208p|j reactions. In the ™''^^Ni-|-^°^Pb reactions the values of ain are close to the 
experimental value of (Ti„ = S-dtfl pb of the ^^Ni+^°^Pb reaction. In the ™''^^Ni+^°^Pb reac- 
tions the cross sections a2n are about 4 times smaller than ain due to smaller values of cTc and 
W,ur- In spite of the large values of Wsur in the reactions ^^'^^Ge+^ospb, 86,88,90,92gg^208pb 
and 88'90'92j^]^_)_208p]^^ ^Yie corresponding values of aER are expected to be smaller than 0.1 
pb due to the very small values of Pcn (Table 2 and Fig. 3). 

In spite of the expected relatively small yields for neutron-rich superheavies the larger 
lifetime of these nuclei will allow a detailed study of their properties. The lifetime of 
molecular-type configurations with an initial DNS in the entrance channel can be studied 
with beams of radioactive nuclei. In reactions with neutron-rich nuclei, a neutron emission 
can occur out of the DNS besides a possible quasi-fission because the characteristic emis- 
sion time becomes comparable with the fusion time. This process decreases the excitation 
energy of the DNS and the fusion probability. The effect of neutron emission from the DNS 
is expected to be important for energies larger than the energy in the 3n channel. With a 
neutron emission from the DNS the fusion process is more complex and has to be studied. 

Intensive beams of neutron-rich nuclei are very useful for producing heavy actinides, e.g. 
Fm as listed in Table 2. In the Pb-based reactions the use of neutron-rich projectiles leads 
to values of aER comparable with evaporation residue cross sections for reactions with the 
stable projectiles. More asymmetric reactions with radioactive beams could be more useful 
in the production of superheavies. 



We thank Profs. D.Habs, R.V.Jolos, G.Miinzenberg, V.V.Volkov, Drs. E.A.Cherepanov, 
S.P.Ivanova and A.K.Nasirov for fruitful discussions and suggestions. G.G.A. is grateful to 
Alexander von Humboldt-Stiftung (Bonn) for support. This work was supported in part by 
DFG and RFBR. 



6 



REFERENCES 



[1] S.Hofmann, Rep. Prog. Phys. 61 (1998) 636. 

[2] G.Miinzenberg, Phil. Trans. R. Soc. Lond. A 356 (1998) 2083. 

[3] Yu.Ts.Ogancssian ct al, Eur. Phys. J A 5 (1999) 63; Preprint JINR, E7-99-53 (1998); 

A.V.Yeremin, V.K.Utyonkov and Yu.Ts.Oganessian m AIP Conf. Proc. 425 (1998) 16. 
[4] D.Habs et ah, NucL Phys. A 616 (1997) 39c. 

[5] G.Miinzenberg et aL, in Superheavy elements. Munich Accelerator for Fission Frag- 
ments. Physics Case and Technical Description, eds. D.Habs et al. (Munich, 1998) 5. 

[6] K.H.Schmidt and W.Morawek, Rep. Prog. Phys. 54 (1991) 949. 

[7] A.Sandulescu, R.K.Gupta, W.Scheid and W.Greiner, Z. Phys. A 283 (1977) 217. 

[8] W.J. Swiatecki, Phys. Scripta 24 (1981) 113. 

[9] P. Frobrich, Phys. Rep. 116 (1984) 337. 
[10] G.Royer and B.Remaud, Nucl. Phys. A 444 (1985) 477. 
[11] C.E.Aguiar, V.C.Barbosa and R.Donangelo, Nucl. Phys. A 517 (1990) 205. 
[12] Y.Aritomo, T.Wada, M.Ohta and Y.Abe, Phys. Rev. C 59 (1999) 796. 
[13] G.G.Adamian, N.V.Antonenko, S.P.Ivanova and W.Scheid, Nucl. Phys. A 646 (1999) 
29. 

[14] V.V.Volkov, in Proc. Int. School-Seminar on Heavy Ion Physics, Dubna, 1986 (JINR, 

Dubna, 1987) p.528; Izv. AN SSSR ser. fiz. 50 (1986) 1879. 
[15] N.V.Antonenko, E.A.Cherepanov, A.K.Nasirov, V.B.Permjakov and V.V.Volkov, Phys. 

Lett. B 319 (1993) 425; Phys. Rev. C 51 (1995) 2635. 
[16] G.G.Adamian, N.V.Antonenko and W.Scheid, Nucl. Phys. A 618 (1997) 176; 

G.G.Adamian, N.V.Antonenko, W.Scheid and V.V.Volkov, Nucl. Phys. A 627 (1997) 

361. 

[17] G.G.Adamian, N.V.Antonenko, W.Scheid and V.V.Volkov, Nucl. Phys. A 633 (1998) 

409; Nuovo Cimento A 110 (1997) 1143. 
[18] R.V. Jobs, A.K. Nasirov and A.I. Muminov, Eur. Phys. J. A 4 (1999) 245; 
[19] E.A.Gherepanov, Preprint JINR, E7-99-27 (1999). 

[20] G.G.Adamian, N.V.Antonenko and Yu.M.Tchuvil'sky, Phys. Lett. B 451 (1999) 289. 

[21] Yu.F.Smirnov and Yu.M.Tchuvil'sky, Phys. Lett. B 134 (1984) 25; O.F.Nemetz, 
V.G.Neudatchin, A.T.Rudchik, Yu.F.Smirnov and Yu.M. Tchuvil'sky, Nucleonic clus- 
ters in nuclei and the many-nucleon transfer reactions (Kiev: Naukova Dumka, 1988). 

[22] A.Diaz-Torres, N.V.Antonenko and W.Scheid, Nucl. Phys. A 652 (1999) 61. 

[23] J.Toke ct al., Nucl. Phys. A 440 (1985) 327. 

[24] D.J.Hinde, et al., Phys. Rev. C 53 (1996) 1290. 

[25] P.Reiter et al., Phys. Rev. Lett. 82 (1999) 509. 

[26] E.A.Cherepanov, in Proc. Int. Symp. on In Beam Nuclear Spectropscopy, Debrecen, 

1984, p.499; E.A. Cherepanov and A.S. Iljinov, Nucleonika 25 (1980) 611. 
[27] P.Moller and R.Nix, At. Data Nucl. Data Tables 39 (1988) 213 
[28] S.Raman et al.. At. Data Nucl. Data Tables 36 (1987) 1. 
[29] A.B. Quint et al., Z. Phys. A 346 (1993) 119. 

[30] C.-C. Sahm et al., Z. Phys. A 319 (1984) 113; Nucl. Phys. A 441 (1985) 316. 

[31] Yu.Ts.Oganessian et al., FLNR Scientific Report 1995-1996 (JINR, Dubna, 1997) 62. 

[32] D.Berdichevsky et al., Nucl. Phys. A 499 (1989) 609. 

[33] H.Gaggeler et al., Z. Phys. A 316 (1984) 291. 



7 



[34] P.Moller et al., At. Data Nucl. Data Tables 59 (1995) 185. 

[35] S.Cwiok and A.Sobiczewski, Z. Phys. A 342 (1992) 203; R.Smolanczuk, J.Skalski and 
A.Sobiczewski, Phys. Rev. C 52 (1995) 1871; R.Smolanczuk, Phys. Rev. C 59 (1999) 
2634. 



8 



TABLES 



TABLE L Calculated maximal energy excess AE^^ above the entrance Coulomb barrier in the 
DNS model and the surplus of energy AB^^p above the Bass barrier extracted from experimental 



data |32|. The number of valence neutron particles or holes to the nearest closed shell in the 
projectile or target is denoted by AA^. 
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TABLE II. Excitation energy E(^^ of compound nucleus, fusion probability Pcn, capture 
cross section ac, surviving probability Wsur, and theoretical <t^^ and experimental evapora- 
tion residue cross sections for reactions leading to the Fm nucleus and Pb-based reactions. The 
experimental data are taken from Refs. |33| and (cold fusion). 
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Fig.1 Adamian et al. 



FIG. 1. Fusion (Pcn) and surviving (Wsw) probabilities as functions of mass number A of 
the projectile in reactions "^Zn+^^^Yb at bombarding energies supplying the maximum yield of 
evaporation residues. 
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Fig.2 G.G.Adamian et al. 



FIG. 2. a) Calculated fusion probabilities Pcn for cold fusion in (HI, In) reactions for the 
projectiles indicated (open circles). The experimental data [||,^ are shown by solid diamonds. For 
the compound nuclei with Z=104-112, the calculations were performed with Q- values from Ref. j27| . 
b) Optimal excitation energies of the compound nuclei. For the nuclei with Z= 113, 114, 116 and 
118, the lower and upper limits of bars were calculated with Q-values from Ref. |27| and Ref. |34|, 
respectively. 
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Fig. 3 Adamian et al. 



FIG. 3. The same as in Fig. 1 for the reactions ^Ni+^°^Pb (sohd lines) and ^Ge+^°^Pb (dashed 
hnes) . 
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